Uranus's anomalously low excess heat constrains strongly interacting dark matter 
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Dark matter annihilations can generate significant amounts of internal heat inside planets if dark 
matter consists mainly of particles with nuclear cross sections in the micro-barn range or larger 
(SIMPs). By considering a detailed model of Uranus's interior, we calculate upper limits on the 
S-wave annihilation cross section for these particles as a function of their mass. These upper limits, 
together with other experimental and theoretical constraints, rule out SIMPs with masses between 
150 MeV and 10 4 GeV. 

PACS numbers: 95.35.+d 



I. INTRODUCTION 



Besides weakly interacting massive particles, many 
other candidates have been proposed for dark matter. 
See pj for a comprehensive review. Strongly Interact- 
ing Massive Particles (SIMPs) have been proposed as a 
solution for the ultra high energy cosmic ray problem 
@ and for the absence of cooling flows 0. Constraints 
from experiments [4|, |5| and from Big Bang nucleosyn- 
thesis and gamma ray observations , limit SIMPs with 
cross sections on protons in the micro-barn range to be 
lighter than 1 GeV or to be heavier than 10 4 GeV. See 
[l4j for accurate calculations of the allowed window below 
4 GeV. We note that, despite the Lee- Weinberg bound 
[l2| . dark matter can consist of particles lighter than 1 
GeV dill. 

In this article we show that annihilations of SIMPs can 
generate significant amounts of heat inside planets. Ob- 
viously, the heat production due to annihilations inside 
the planet has to be smaller than the total internal heat 
production. This provides us with a new way to con- 
strain SIMPs. Large planets with very low excess heat 
pose the strongest constraints. Uranus's internal heat 
production is atypically small, only about a tenth of the 
similar sized planet Neptune [llj. In this idea was 
used to exclude the H dibaryon model and in H new 
constraints have been obtained for other baryonic dark 
matter candidates. 

The capture rate of SIMPs is negligible unless 
^criipianet >> 1, where n is the density of nucleons in 
the planet, i? p i a nct the radius of the planet and a the av- 
erage elastic cross section. This means that a needs to 
be larger than about 10 -32 cm 2 for SIMPs to generate 
a significant amount of heat inside large planets such as 
Uranus. We can crudely estimate an upper bound on 
the heat production, W mH , by assuming that the entire 
dark matter flux is transferred into heat. Assuming a 
relic dark matter density of 0.3 GeV/cm 3 and a typical 
velocity of 270 km/s of a dark matter particle w.r.t. the 



solar system we get: 

W max « 7ri? 2 lanct x (0.3 GeV cm" 3 ) x (270 km/s) 

« 4 ( S^l) x 10 15 W 
V 10 4 km / 



(1) 



In case of Uranus, i? p i an ct = 2.509 x 10 4 km and thus 
Wmax ~ 2.5 x 10 16 W, which is much larger than the 
observed heat production of (3.4 ± 3.8) x 10 14 W [TJ. 

To calculate the actual heat generated by dark matter 
annihilations in a SIMP dark matter scenario, we have 
to take into account that, because of the mass mismatch 
with nuclei, SIMPs need many collisions to lose enough 
energy to become captured inside a planet. For low mass 
SIMPs this means that a significant part of the incident 
flux of SIMPs will end up being reflected back into space. 
Also, if the annihilation cross section is low enough, most 
SIMPs will eventually evaporate rather than annihilate. 
The evaporation rate depends not only on the mass of 
the SIMPs but also on the elastic cross section. Particles 
escaping from the gravity of a planet typically escape 
from one scattering length below the surface. Particles 
with lower elastic cross sections will thus escape from 
deeper below the surface, where it is hotter. This means 
that these particles will have higher escape rates. 

We will use a detailed model of Uranus's interior to 
take these effects into account for SIMPs lighter than 1 
GeV interacting with spin independent interactions with 
nuclei. Upper limits on the S-wave annihilation cross 
section will be derived as a function of the mass and the 
elastic cross section on protons. 



II. CAPTURE RATE OF LIGHT SIMPS 

The velocity distribution, / (v), of halo dark matter 
particles is usually assumed to be Maxwellian 0: 



f(v) 



exp 



(2) 
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Here vq ~ 220 km/s 0| is the local circular speed. This 
velocity distribution has to be cut off at the escape speed 
of the galaxy. Since particles in the tail of the velocity 
distribution will make only a small contribution to the 
captured flux, we will ignore this. The solar system is 
moving around the center of the galaxy with a speed 
of about vq. This means that the velocity distribution 
w.r.t. the solar system is anisotropic and rather compli- 
cated. However, since we have to integrate the flux of 
the SIMPs over the surface of a planet, we can replace 
the actual velocity distribution by the one averaged over 
all directions. This averaged velocity distribution, / (y), 
is: 
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A SIMP entering a planet will typically undergo multi- 
ple scatterings before it is either reflected back into space 
or is thermalized and captured. In [h} this problem is 
investigated for the Earth for dark matter particles heav- 
ier than 2 GeV using simulations. In the Appendix we 
fit these results to an approximate expression obtained 
by treating this problem in the continuum limit, where 
it becomes a diffusion problem. We find that the capture 
rate, denoted as lZ ca , P t, is: 



n, 



capt 



4nR 



planet n simp 



R (v) nf («) v 3 dv (4) 



where n s i mp w 0.3 GeV/(m s i mp cm 3 ) is the halo number 
density of SIMPs, v csc is the escape velocity from the 
planet and R (v) is the probability that a SIMP moving 
into the planet with speed v is captured: 
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1.27 x 
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Here m c g is an effective nucleon mass, defined as: 



(m e ff)~ 



(5) 



(6) 



where rrij denotes the mass of a nucleus of type j and 
fj is the relative probability that a SIMP will scatter off 
a nucleus of this type, fj is thus proportional to rij(Tj, 
where rij is the number density of nuclei of type j and 
<jj is the elastic cross section with nuclei of that type. 
For spin independent interactions the Gj can be related 
to the elastic cross section with protons as follows ^{| : 



Mi 



42^ 

-rij °prot ■ 



(7) 



Aj is the 

number of nucleons in nuclei of type j, fj,j is the relative 



Here <7 prot is the cross section with protons, ^ 



mass of a SIMP and a nucleus of type j and /x prot is the 
relative mass of a SIMP and a proton. We will assume 
that SIMPs undergo most of their scatterings before be- 
ing captured in the upper 20% of Uranus's atmosphere. 
There 93 % of the nucleons are protons and about 7% 
are helium nuclei. According to JBJ, the effective nuclear 
mass is: 



0.93 1 



0.93 1 



+ + 1.12 ( 1 + ^isp y 

( 1 + ^h^) 2 + 0.28 f 1 + 



(8) 



And the average cross section, cr, is: 



0.93+1.12 



4m p 



^prot • 



(9) 



III. 



EVAPORATION AND ANNIHILATION OF 
SIMPS 



To estimate the evaporation and annihilation rates of 
captured SIMPs, we need to know the density profile 
of the SIMPs. As long as temperatures don't change 
too rapidly over the course of the length of one mean 
free path of a SIMP, we can assume that the SIMPs are 
in local thermal equilibrium with the nuclei in Uranus's 
interior. Local thermal equilibrium will break down high 
in Uranus's atmosphere where densities are low. Since 
only a small fraction of the SIMPs is expected to reside 
there, this can be neglected in a first approximation. 

The density profile, n s i mp (r), follows from demanding 
that the SIMPs are in hydrostatic equilibrium: 



ap (r) 



Ttsimp (0) T (0) 

T(r) 



exp 



>G 



M(y) 
y 2 T (y) 



dy 



(10) 

Here G is the gravitational constant and M (y) is the cu- 
mulative mass of Uranus up to radius y. Even though it is 
impossible to measure the density and temperature pro- 
file of Uranus, reasonable accurate models can be made 
based on reasonable assumptions about the interior com- 
position, the measured heat flux, the total mass and grav- 
itational moments |16|. The density profile and temper- 
ature profiles based on such a model are shown in Fig. 
and Fig. [2 respectively [l^. In Fig. [3| the density profile 
of 0.5 GeV SIMPs is shown. 

The evaporation rate, lZ cvap , can be estimated by cal- 
culating the upward flux of SIMPs from one scattering 
length below the surface |20j: 



(11) 
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FIG. 2: Interior temperature (K) as a function of the radius 
in meters. 



Here is: 



/ n (r) adr = 1 



(12) 



and v osc (r* ) is the escape velocity at r = r, . 

The annihilation rate, 7?. a , m , can be expressed as: 



CI XTlgiYap 

6k 



a J d xn simp (x) + ——^b j d xn simp (x) T (x) 

(13) 

Where we have assumed that cr ann is well approximated 
by the S-wave and P-wave contributions: 



b—. 



(14) 



IV. CONSTRAINTS ON cr a nn 

Besides satisfying the constraints from Uranus's excess 
heat, the annihilation cross section has to be consistent 
with the SIMPs being thermal relics. According to 01 
for particles in the MeV-GeV range: 



10" 27 cm 3 



< 



< 



1CT 26 cm 3 



/s (15) 



during the freeze-out epoch. For m s i mp in the range of a 
few hundred MeV to 1 GeV the freeze-out temperature 



is about 



'/ (20k). This means that at freeze-out 



farm ~ a + 6/3. The lower limit for the annihilation 
cross section during freeze-out thus implies that: 



+ b/3> 10 



-27 



Vs. 



(16) 



Uranus's excess heat, as measured by Vo yag er 2 dur- 
ing its 1986 fly-by, is (3.4 ± 3.8) x 10 14 W [llf. Models 
of Uranus's atmosphere based on the observed temper- 
ature profile in Uranus's upper atmosphere suggest that 
Uranus's internal heat is close to Voyager's 1-cr upper 
bound 17]. We shall consider SIMP models that predict 
that more than 10 15 W of energy is generated by anni- 
hilations as excluded. Although only part of the energy 
generated by annihilations can contribute to the internal 
heat, it is unlikely that a significant fraction of Uranus's 
internal heat is produced by exotic processes such as an- 
nihilations of SIMPs. 

To estimate the energy produced by SIMP annihila- 
tions, consider the time evolution of the total number of 
captured SIMPs, Af (t): 



dM{t) 
dt 



= n 



capt 



cvap 



(17) 



It follows from IjlUfl . (fill) and (|13fl that there exists func- 
tions a (m s i mp , (Tprot) and /3(m s i mp ) such that: 



n, 



cvap 

V 

/v ann 



a (m simp , cr prot ) A/", 

^(TOsimp)^ 2 . 



(18) 

(19) 
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Note that j3 (m s i mp ) is a linear combination of the pa- 
rameters a and b. We set t — at the time of formation 
of the solar system, about 4.5 billion years ago. Thus 
TV (0) = and (|T7|> implies that: 



Af(t) 




capt tanh ( ^-Ja 2 + 4/3ft capt (t + C) 



(20) 



where 



C 



^/a 2 + 4/3ft capt 



: arctanh 



Q 



4(311, 



capt 



(21) 

The total energy generated by SIMP annihilations now 
follows by substituting this in l|19l) and putting t = 4.5 
billion years. 

Demanding that m simp c 2 [ (iN ' (t) 2 < 10 15 W yields the 
following constraints on the annihilation cross section. 
For every m s i mp there exists a critical value, er cr it for 
Oprot- If c prot is smaller than er cr it, there are no con- 
straints on a and b except for the usual thermal relic con- 
straints. If cr pro t is larger than cr cr it, a has to be smaller 
than 10 -27 cm 3 /s. According to l|16fl this means that 
b can no longer be zero. The upper limit for a quickly 
reaches a terminal value, a max , for cr pro t larger than cr cr it. 
To a good approximation this happens for <7 pro t > 3cx, 



^crit ■ 



The values for a rr \t and 



for 



''Simp 



GeV to 0.3 GeV are shown in Tableffl 



in the range 0. 





fcrit 


timax 


GeV 


lC)- al cm 2 


10- a4 cm a /s 


0.8 


1.6 


9.0 


0.7 


2.9 


7.0 


0.6 


5.6 


8.8 


0.5 


12 


11 


0.4 


28 


15 


0.3 


86 


17 



TABLE I: Constraints on the annihilation cross section as a 
function of m s i mp and cr prot . If a pTOt > <7 C rit, a has to be 
smaller than 10 -27 cm 3 /s and thus b can no longer be zero. 
If Cprot > 3a cr it, a has to be smaller than a ma x. 



Since it is to be expected that the S-wave cross section 
should not differ from the elastic cross section by more 
than a few orders of magnitude, these results rule out 
SIMPs with cross sections larger than er cr it. 

For m simp < 0.3 GeV, the density of SIMPS becomes 
large at the outer regions of Uranus's atmosphere. To de- 
termine cr cr it and a max with reasonable accuracy would 
require using detailed models of Uranus's upper atmo- 
sphere. We estimate that for m s i mp < 0.15 GeV, SIMPs 



cannot be constrained and that at m smip = 0.15 GeV, 
CT pro t has to be smaller than 10~ 27 cm 2 . For m simp > 0.8 
GeV, cr cr it becomes so low that local thermal equilibrium 
starts to break down. It is clear, however, that SIMPs in 
the mass range 0.8 - 1 GeV are excluded as well. Above 1 
GeV known experimental constraints exclude SIMPs up 
to masses of about 10 4 GeV 



V. CONCLUSIONS AND OUTLOOK 

We have shown that Uranus's low excess heat rules 
out SIMPs as dark matter candidates for 150 MeV < 
TOsimp < 1 GeV. SIMPs with masses in the range 1-10 4 
GeV were already ruled out ||. It would be interesting 
to constrain SIMPs above this mass range. Preliminary 
results indicate that for 10 4 GeV < m sirop < 10 10 GeV 
the constraints on the elastic cross section on protons 
posed by Uranus's excess heat are stronger than the re- 
sults of Limits on baryonic dark matter candidates 
0,1- 



have been obtained in 
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APPENDIX A: CALCULATION OF THE 
CAPTURE RATE 

In this appendix we will derive Eq. 0] for the capture 
rate of light SIMPs. We will first derive an approximate 
analytical expression and then modify this somewhat to 
fit the results of simulations given in ^4[. The analyt- 
ical expression will be derived by estimating how many 
collisions are typically needed for a SIMP to slow down 
below escape speed. We then evaluate the probability 
that a SIMP entering the planet stays below the surface 
after this number of collisions. 

We will assume that the mass of the SIMP, m s i mp , is 
much less than the typical nucleon mass and that elastic 
cross sections with nuclei are isotropic. Under these con- 
ditions a SIMP inside a planet will perform an isotropic 
random walk. After each collision the fractional energy 
loss is a random number, uniformly distributed in the 



5 



interval /: 



0,1 



1-r 
1+r 



(Al) 



Here r 



j/wnudeon- This is true for general r, not 



just when r << 1. We can thus easily estimate how 
many scatterings are typically necessary to reduce the 
energy from ^m simp w 2 to |m S i mp (v csc ) , where v, ' 



is 



the escape velocity from the planet. Denoting the ratio 
of the energies after and before the j-th scattering by 7^ , 
we can estimate the energy after N scatterings, En, as 
follows. We first write: 



JV 



E N = E Q JJ jj = E Q exp 



3 = 1 



N 



£ lo s fa) 



(A2) 



Let's denote the relative probability that the SIMP will 
scatter off a nucleus of type p by f p . Typically, a fraction 
f p of the 7j in l|A2() is an energy ratio before and after 
scatterings off nuclei of type p. These 7's are uniformly 



distributed in the interval 



V l + r P ) 



1 



where r p is the 



ratio of the mass of the SIMP and a nucleus of type p. 
The average of log (7) is thus: 



(log (7)) = 1 + 



(1 



?>) 2 log 



2r„ 



(A3) 



By replacing the log (7) 's by their averages in (|A2|> the 
following approximation is obtained: 




En ps Eq exp 



It follows from (|A4(I that the number of collisions needed 
for a SIMP to be slowed down from a speed of v to the 
escape velocity v csc , N (v), is typically: 



N(v) 



log 



fc simp 



^csc 



log 



l+r p 



(A5) 



Here we have used that r p << 1 and defined 



-1 — \ 



(A6) 



where m p is the mass of a nucleus of type p. 

Next we will evaluate the probability that a SIMP at 
some distance below the surface will remain below the 
surface after N (v) scatterings. As long as onR >> 1, we 
can treat this problem as a random walk in a half infinite 
space with absorbing boundary conditions. If N (v) is 



not too small, the random walk can be treated in the 
continuum limit where it becomes a diffusion problem. 
To find the continuum limit in this case, consider the 
evolution of the probability distribution of a SIMP in an 
infinite medium. The probability density Q (x, xl) that a 
SIMP will scatter at position x if the previous scattering 
was at position xl is: 



! (x, xl) 



an exp (— no \x — xf\) 
4-7T (x — xl) 2 



(A7) 



The probability distribution for the SIMP after N scat- 
terings, Pjv (x) , evolves according to: 



(A8) 



P N+ i (x) = / d 6 xiP N (xi) Q (x, xi) . 



It follows from (| AT|> and l|A8(l that the Fourier transform 
of Pn (x), defined as: 



p N (k) = I d 3 xP N (x) exp (-27rifcx) 

satisfies 

P N+1 (k) = A (k) 



_27r|fc| 

For large N this becomes: 



no' j 2tt \k\ 

arctan 
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1 N 



N /2nk 
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P N+1 (k) « Pi (k) exp 
This implies: 

(x) w / (fxiPi (xi) Q N (x, xi) 



Pi 



N+l 



(A9) 



(A10) 



(All) 



(A12) 



where Qn (x, xl) is: 



)jv (x, xl) 



3n 2 o 2 
AttN 



exp 



3n 2 o 2 (x - x/y 
AN 



■;Ai3i 

To find the evolution of the probability distribution of 
a SIMP in a half infinite space, we must put Pn (x) — at 
the boundary. For problems that can be formulated en- 
tirely in the continuum limit this can be exactly achieved 
using the method of images, see e.g. |l8|. According to 
this method, if g (x) is the initial probability distribution 
of the SIMP defined in the half infinite space, we should 
substitute for P x (x) in (|A12|) : 



P 1 (xi)=g(xi)-g(R(xi)) 



(A14) 



where R (xl) represents the vector obtained by reflecting 
xl in the boundary of the half space. We expect that this 
method will yield a reasonable approximation. To obtain 
the number of SIMPs captured by a planet per unit time 
and per unit area we need to substitute for g (xl) the 
number density of SIMPs scattering for the first time at 



6 



xi per unit time and per unit volume in velocity space, 
01 : 

<t>! (a,v) = n simp exp (- ™^ ^ naf(v) \v\ . (A15) 

Here a is the depth below the surface, n s i mp the halo 
number density of SIMPs and 6 is the angle of the velocity 
w.r.t. the normal of the surface. The captured flux of 
SIMPs, ^" captur ed; can now be obtained by substituting 
(|A15(I in l|A14(l and the resulting expression for P\ in 
(|A12|I . By integrating this over xi and x over the half 
infinite space below the surface we obtain to lowest order 
in l/N{v): 

(j)(v)^ ^=n simp N{v)'Kf{v)v 3 . (A16) 

V07T 

Here 4> (v) dv is the flux of captured particles with original 
speeds between v and v + dv. For speeds v ~ u CS c the 
above derivation is not very accurate because then N (v) 
is not large. However, since / (v) is small in this case, 
this isn't a problem. We shall neglect the flux of particles 
moving slower than v csc . Integrating (|A16(1 and using Eq. 
IA5l for N (v), we obtain for the flux of captured particles: 

/>oo 

•^captured = "simp / R (v) TT f (v) V 3 dv . (A17) 



Here 

R(v)~*[^L g (^)Y L2 (A18) 

V37T V m cff \ V V csc J J 

is the probability that a particle moving with speed v 
into the planet is captured. 

In the capture of strongly interacting dark mat- 
ter particles heavier than 2 GeV by the Earth is investi- 
gated using simulations. We find that Eg. IA18I predicts 
a captured flux that is about a factor of 2 too low. The 
likely source of this difference is the improper use of the 
method of images to impose absorbing boundary condi- 
tions. Changing the prefactor of — 1= to 1.27 gives an 
almost perfect fit to the captured fraction of the flux in 
the mass range 2 GeV-10 GeV. For Earth, m e ff ~ 21.5 
GeV, about ten times as high as for Uranus. When rn s i mp 
is below 10 GeV, the ratio m s i mp /rn e ff for Earth is thus 
in roughly the same range as the ratio m smip /77i e ff for 
Uranus when m s i mp is below 1 GeV. We thus expect that 

i2( w )«1.27 [™^Lg(^-)) 5 (A19) 

V m cff \ \V CS c/j 

is a reasonable approximation for the capture probability. 
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